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Abstract 

We study the 2- loop resummed propagator in hot <? 2 </> 4 theory. The propa- 
gator has a cut along the whole real axis in the complex energy plane, but 
for small g, the spectral density is sharply peaked around the plasmon. The 
dispersion relation and the width of the plasmon are calculated at zero and 
finite momentum. At large momenta the spectral width vanishes, and the 

plasmon looses its collectivity and behaves like a non-interacting free particle. 
PACS numbers: 12.38Cy, 12.38Bx, ll.lOWx, 12.38Mh 
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I. INTRODUCTION 



At high temperature T and low momenta p < T the effective degrees of freedom in any 
field theory are collective modes, bosonic "plasmons" |T|,0] and fermionic "plasminos" 0,ffl]. 
Even if the fundamental field theoretical degrees of freedom were massless, the collective 
modes possess a finite "thermal" mass which is generated dynamically by the interactions 
among the fundamental degrees of freedom. The thermal masses regulate (some of) the 
severe infrared divergences which the massless theory would otherwise develop at finite 
temperature due to the singularity of the Bose distribution at zero momentum |5]U. The 
strong rise of the Bose distribution at small momenta causes, for example, the transport 
properties (viscosity, heat conduction) of a hot plasma of massless fields to be dominated by 
the interaction between low-momentum collective modes The determination of their 

dispersion relation and their collisional width is therefore a necessary prerequisite for any 
microscopic calculation of the transport coefficients, for example using the Kubo formulae 

It is known that a consistent determination of the plasmon dispersion relation and width 
in a massless field theory requires the resummation of certain loop contributions to the 
propagator, the so-called "hard thermal loops" (HTL) ||. In massless gauge theories, this 
resummation becomes non-trivial beyond leading order, due to the complicated analytical 
structure of the HTLs (i.e. their momentum dependence and logarithmic branch cuts from 
Landau damping |L]|| ) . In massless scalar 4 theory, the resummation is much easier 0,^] , 



since the leading HTL (the tadpole diagram) is just a momentum-independent real constant, 
and it can therefore be relatively easily be carried beyond leading order ||11|| . Still, when 
we started this work, no calculation of the collisional width of finite momentum plasmons 
existed. Since the latter is needed, however, for the calculation of transport coefficients 
from the Kubo formulae (previous calculations H3|,|l"3f used, for lack of better knowledge but 



without justification, the zero momentum limit of the plasmon width), we present in this 
paper a calculation of the resummed scalar field propagator in g 2 (p 4 theory at the 2-loop 
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level, including real and imaginary parts in 4-dimensional momentum space. While this 
work was in progress, Jeon [[12] presented a calculation of the plasmon width at zero and 
finite momentum on the plasmon mass shell. Our results agree with his and extend them 
into the off-shell domain. 

This paper is organized as follows. In Section [TT] we will give a brief review of results 

we 



from Ref. |TT[ on resummation in hot theory which we will need later. In Section [Ti I 
present the calculation of the imaginary part of the 2-loop self-energy and investigate the 
analytical structure of the full propagator at two-loop order. The properties of plasmons at 



rest and at finite momentum are studied analytically and numerically in Sections [III A| and 



1 1 [B|, respectively. Our conclusions are summarized in Section |V[ Appendix |A] contains 



some technical steps of the calculations in Section LI I 



II. RESUMMATION AT HIGH TEMPERATURE 

2 

For a hot scalar field, we consider the following Lagrangian with A = < 1: 

£o = \{d»<P? + 9 ^. (2.1) 
At the tree level it describes massless scalar fields with a quartic self interaction. The induced 



thermal mass resulting from the single HTL (tadpole) is of order g and reads |ITj| . |TT 

™ t 2 h = ^j- • (2-2) 

The effects from this thermal mass can be resummed by defining an effective Lagrangian 
through 

£ = (Aj + \ m \^ - \ m \^ 2 = £eff ~ 2 TO th0 2 > (2-3) 

and treating the last term as an additional interaction. This effective Lagrangian defines an 
effective propagator with thermal mass m t h 



where in the imaginary time formalism K 2 = k$ + k 2 , ko = 2nirT, n = 0, ±1, ±2, .... It 
represents the infinite sum of iterated HTLs. Since m t h leads to a qualitative modification of 
the resummed effective propagator for soft momenta, p < gT, its inclusion via resummation 
is necessary to avoid infrared divergences at finite temperature ||. Vertex corrections from 
HTLs are down by a factor g relative to the bare vertex g 2 and don't qualitatively change 
the latter. In a perturbative expansion of 4 theory we can thus keep using the bare vertex 
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One can now recalculate the 1-loop self-energy using this resummed effective propagator. 
This amounts to computing the contribution of all "cactus diagrams" or "superdaisies" in the 



original Lagrangian [10] each of which is infrared divergent. Using dimensional regularization 



in the MS renormalization scheme, the diagram in Fig. |T]a gives 

)+0(g 5 \ng), (2.5) 



V - 9 2 T 2 ( , 3?77,th 3 fm th \ 2 



ln (^) +7 * 



where /i is the scale parameter in dimensional regularization and 7^ = 0.577216... is Euler's 
constant. The resummation of an infinite series of infrared diagrams reflects itself through 
a g 3 correction to the leading term — g 2 T 2 /24 = —m 2 h ; it is non- analytic in the coupling 
constant A = g 2 /24 and for g <C 1 dominates any finite genuinely perturbative corrections 
of order g 4 . In this sense resummation is essential. 

Including also the contribution from the new 2-point interaction in ( |2.3p , Fig. [TJb, the full 
1-loop self-energy can be written as £1 = £' x +m 2 h . It is purely real. Thus, to order g 3 , the 



full propagator has only two poles at po = ±y P 2 + f^l where m p is the 1-loop resummed 
plasmon mass: 



m p = m th J 1 - . (2.6) 



Following Ref. |11|| , we can now improve the resummation procedure by including into the 
effective Lagrangian all 1-loop effects which are independent of the regularization scheme, 
i.e. by writing 

C = (£0 + \my) - (2.7) 



and expanding into effective scalar propagators with an effective plasmon mass m p . Then all 
factors m t h in Eq. ( |2.5| ) are replaced by m p , and from the full 1-loop self energy = T! x -\-mi 
all terms up to order g 3 cancel, leaving a leading contribution of order g 4 which is real. 

At 2-loop order, the Feynman diagrams shown in Fig. |2] contribute. Their leading con- 
tributions are also of order g 4 . Diagram 0c arises from the new 2-point interaction of the 
Lagrangian ( |2.7|) and contributes at the same order. 



The contribution E' 2 of the diagrams ||a and |]c is easily evaluated. It is purely real. As 
a consequence of resummation, all leading g A terms cancel, and the result is found to be 

V 2 = 0(jg 5 ]ng). (2.8) 

A convenient way to calculate diagram |2|b is to use the Saclay method Hl4| . Its contribution 
can be expressed as 

£2(^0, P) = / d[k, g] {S(E k , E g , E r )[(l + / fc )(l + + f r ) - f k f q f r ] 

+ S(E k , E q , -E r )[(l + f k )(l + f q )f r - f k f q (l + f r )} 
+ S(E k , -E g , E r )[(l + f k )f q (l + f r ) - f k (l + f q )f r ] 
+ S(-E k , E g , E r )[f k (l + f q )(l + f r )-(l + f k )f q f r }} , (2.9) 



with the same notation as in Ref. [II], namely 



g A j J i£ d D ~ 1 k d D ~ 1 q 1 

d[k > q] = — Wy^ Wy^ mE-E r ' (2 - 10) 

S(E k ,E ,E r ) = , (2.11) 

1 k) qi r> ip + E k + E q + E r ipo -E k -E g -E r ' 1 ; 



where 



r=|k + q-p|, (2.12) 

E? = l 2 + m 2 p , l = k,q,r, (2.13) 

/, = r , l = k,q,r. (2.14) 

exp(Ei/T) - 1 

D = 4 — 2e is the the number of space-time dimension introduced in dimensional regu- 
larization. Please notice that now the one-loop resummed mass m p occurs in all on-shell 



energies. The real part of £ 2 ' contains an ultraviolet divergence. A detailed discussion of the 
subtraction of this ultraviolet divergence by a counter-term in the Lagrangian and the final 



expression for the real part of the renormalized self-energy can be found in Ref. ||11|| . We 
can see from Eq. ( |2.10| ) that Re £ 2 ' is of the order g A and depends on the mass parameter \x 
of the dimensional regularization. 

The complete 2-loop expression for the real part of the self-energy is given by 

Re£(ip ,p) = £1 + £ 2 + Re£ 2 '(ip ,p) , (2.15) 



with the last term taken from Ref. |[L1|| . A detailed calculation of the imaginary part of the 
self-energy will be given in the following section. 



III. FULL PROPAGATOR AND SPECTRAL FUNCTION OF THE PLASMON 

The leading contribution to the imaginary part of the self energy comes from the 2- 
loop diagram |2|b. Making an analytical continuation ip Q = u + irj, 77 = + , and using 
I /{A + irj) = V{l/A) - m8{A), we obtain from Eq. Q 

Im£(u;,p) = Im £'2(0;, p) = — tt J d[k, q] (3.1) 
{[8(u + E k + E q + E r ) - 8{u -E k -E q - E r )][(l + f k )(l + f q )(l + f r ) - f k f q f r ] 
+[8(oj + E k + E q - E r ) - 8(u -E k -E q + E r )][(l + f k )(l + f q )f r - f k f q (l + f r )] 
+[8(u + E k -E q + E r ) - 8(u -E k + E q - E r )][(l + f k )f q (l + f r ) - f k (l + f q )f r ] 
+[8(oj -E k + E q + E r ) - 8(u + E k -E q - E r )][f k (l + f q )(l + f r ) - (1 + f k )f q f r }}. 

This gives the rates for physical decay and scattering processes IfOifl : The first line corre- 
sponds to the decay of one off-shell into three on-shell <fr mesons and (with opposite sign) the 
inverse process; the two 5-functions represent the decay of particles with negative and pos- 
itive energy, respectively. The remaining lines in Eq. ( pj.l[ ) correspond to Landau damping 
via scattering of the off-shell meson with on-shell particles from the heat bath, forward and 
backward scattering process contributing with opposite signs. Landau damping requires the 
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presence of thermally excited states and vanishes at zero temperature. In each term every 
ingoing meson from the heat bath receives a thermal weight fa and each outgoing one a 
Bose-enhancement factor (1 + fa). 
It is easy to show that 

Re E(— uj, p) = Re £(&>, p) and Im S(-oj, p) = -Im p) , (3-2) 

so we restrict ourselves to positive energy u > 0. After suitably relabeling the integration 
variables, the imaginary part of the self-energy ( |3.1|) for u > can be reduced to 

ImS(w,p) = Imgi(w,p) + lmg 2 (uJ,p) (3.3) 

with the 3-body decay and Landau damping contributions, respectively: 

Im^p) = 7r(e^ T -I) J d[k,q]f k f q f r 5(u -E k -E q - E r ) , (3.4) 
hng 2 (u,p) = 3ir(e^ T -l) J d[k, q](l + f k )f q f r 5{u + E k - E q - E r ) . (3.5) 



The ^-functions in Eqs. ( |3.4j ) and ( |3.5| ) arise from energy conservation and constrain the 
integration over k and q to the kinematically allowed domain. The evaluation of the resulting 
kinematic limits of the integration variables is the most demanding part of the following 
calculation and is shortly discussed in Appendix |A[ Before going into the technical details 
of evaluating the integrals, let us summarize some relevant analytical properties of the 
propagators and self-energies. In the imaginary time formalism, the full propagator 

fl(vo,p)= -p§+P 2 +4-£(* PO ,p) (3 - 6) 

has the spectral representation 



oo 



£>(/>„ P)=/ du'&^- r (3.7) 



oo 



The spectral density satisfies 



p) = \D(u + ir], p) - D(u - ir], p) 

1 ImE(w,p) 



7r [to 2 — p 2 — m p 2 + ReS(c<j, p)] 2 + [Im p)] : 



(3.8) 



It is an odd function of uj, see Eqs. ( |3.2|) : 

p(-w,p) = -p(w,p) or p(w,p) = e{u)p(\u\,p) , (3.9) 

where e(ou) = sgn(u;). Thus the full propagator D(ipo, p) has the following analytical struc- 
ture: If ImS ^ 0, the full propagator has a cut along the real axis in the complex uj plane, 
with a discontinuity across this cut given by the spectral function. If Im £ goes to zero, the 
identity 

- lim - V = 5(A) (3.10) 
7i v-o A 2 + r] 2 K ' 

shows that the spectral density can be expressed as a 5-function at uj 2 = p 2 +m?— Re £(u;, p); 
the corresponding analytical structure of the full propagator D(ipo, p) are poles in the com- 
plex uj plane given by the zeroes of the 5-function. 

From D(ip , p) the full retarded propagator D(lo + ir], p) is obtained by analytical con- 
tinuation ipo — > uj + ir\. The full propagator in the real time formalism can be written as 



15 



G(uj, p) = D(uj + ie(uj) V ) - 27rz g^M^ . (3.11) 



All of these full propagators have the same analytical structure. 
Defining the damping rate as 



7(",P) = (3-12) 



and substituting this into Eq. (ft.8|), we see that if 7(0;, p) <^ J p 2 + m v 2 — Re Y,(uj, p), then 
p(u,p) is strongly peaked at uj 2 = p 2 + m p 2 — ReE(a;,p). The solution of this dispersion 
relation defines a weakly damped single-particle state which is interpreted as a quasiparticle, 
the so-called plasmon. Denoting the solution for uj by e(p) and defining T = 27(e(p),p), 
the spectral density can in this limit be expressed approximately in the form of a relativistic 
Breit-Wigner function with full width T: 

1 r 

p(uj,p)tt-— — — - (3.13) 

71 (uo z — e (p)) + <jJ i 
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Finally, the spectral function fulfills the sum rule [|T6 



roc 

2 duu p(u, p) = 1 
Jo 



(3.14) 



In the following two subsections we now investigate in detail the properties of plasmons 
in 4 theory, both at rest and at finite momentum p relative to the heat bath. 



A. Plasmons at rest, p = 

For p = the evaluation of the imaginary part of the 2-loop self-energy ( |3.3|) can be 
carried out in the following manner: for the q-integration we choose the z-axis along k and 
define t = cos(k, q). The integrals in Eqs. (|3.4|) and (|3~5|) thus take the form 



7r J d[k, q] 



9 I J7„ r„2 / j„„2 / 1 



dkk 2 dqq 2 dt , (3.15) 

48(2tt) 3 7o Jo HH J-i E k E q E r ' v ; 



where E r = ^jk 2 + q 2 + 2kqt + m 2 . The energy conserving (^-functions in Eqs. (|3.4|J3.5|) will 
be used to perform the t-integration by using 

dE r kq 

■w = i- (316) 

The condition that the zeroes of the 5-function fall into the interval — 1 < t < 1, i.e. 



(k — q) 2 + m p 2 < E r < yj \k + q) 2 + m p 2 , leads to restrictions of the integration regions in 
k and q. These restrictions are worked out in Appendix [A|. For the 3-body decay contribution 
we find 

g'T 2 



Im^,0) = ^-3m p )|-— (W T -l) 



384tt 3 

x/ dvf(v)f(w-v) In r , 3.17 

•/« V/( e (^))/(^(yi) +v- w)J 

where 

a = — , u> = -, a=~, y=^, e(z) = Vz 2 + a 2 , /(v) = -^—j . (3.18) 

The specific values x\ and yi 2 are given by Eqs. ( |A7|) and ( |A4|JA5| ) for p = 0, respectively, 
divided T. 



Similarly, the Landau damping contribution Im g2(u,0) is given by 
9 4 T 2 



lmg 2 (uj,0) 



1287T 3 



(e"' T - 1 



8(m P 



x(l + f(v))f(w + v) In 



-co) dv + 6{uj — m p ) 

Je(xl) 

7 (£(2/3)) / (£(2/4) -«-«>)' 



(3.19) 



J (e(y 4 )) f (e(y 3 ) - v - w) J ' 

with £4 and 2/3,4 given by Eqs. ( |A17| ) and ( |A12j , |A13] ) for p = 0, respectively, divided by T. 

While Im<7i(u;,0) vanishes below the 3-particle threshold Eq = 3m p (see Eq. ( |A2|) ), the 
Landau damping contribution Im g^O^O) (which arises only at T 7^ 0) is nonzero for all 
positive values of lu. Due to the symmetry (|3.2|) , at 2-loop order the full propagator thus 
has a cut along the whole real lu axis. This is qualitatively different from the simple double 
pole structure at 1-loop order. 

On the plasmon mass shell lu = m p only Landau damping contributes to the imaginary 
part. From Appendix [X] we find that in this case B(m p , k) = k 2 and thus ys = 0, = x as 
well as x\ = 0. Substituting this into Eq. (|3.19| ) and using Eqs. ( |3.3| ) and ( |3.12j) , we obtain 



the on-shell plasmon damping rate (see also |I2]) 



7(m p ,0) 



ImS(m p ,0) g A T 2 



2m r 



256n 3 m, 



-L 2 (e- a ) 



9 3 T 



where L 2 (z) is the Spence function defined as 



64V247T 



[l + 0(glng)} 



(3.20) 



Lo(z) = - 



Z dt l -^1 



(3.21) 



The leading term in the last expression of ( 3.20|) was also found in Ref. [ I I 



Im E is finite and does not require any regularization. For dimensional grounds the width 



7(m p , 0) for plasmons at rest is thus proportional to the temperature T. From ( 3.2C ) we 
further see that for small g the ratio j(m p , 0)/T is proportional to g 3 . 

Using Eq. ( |3.8| ) we can also evaluate the spectral function. We notice that, as discussed 
in Section [TT], the 0(g 4 ) term in ReE contains the renormalization scale \x which should 
be absorbed into the renormalized coupling constant. But since the 0(g 4 ) corrections have 
only a minimal effect on the position of the peak in the spectral function, and all 0(g 3 ) 
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terms have already been included in the plasmon mass, we can to very good approximation 
neglect ReS. In this approximation the plasma frequency 



. 1J = £(p = 0) = rop = ^l--^L (3.22) 

is strictly proportional to T. An additional logarithmic temperature dependence through 
the running of the coupling constant can only enter in higher order. 

In the following we comment on some further properties of the plasmon based on numer- 
ical results. In Fig. |]a we show for m^/T = 0.1 the 2-loop spectral function together with 
a line indicating the pole position at one-loop order (including HTL resummation). Note 
that there is no visible shift of the peak due to the higher order corrections; the shift of 
the peak below the value of m th is entirely due to the 1-HTL- resummation, see Eq. ( |3.22|) . 
The major new feature at 2-loop order is just the finite width of the plasmon. Although 
the 2-loop spectral density is strictly positive for all values of u, it approaches zero very fast 
as one moves away from the plasmon peak. Only in the neighborhood of the 3-body decay 
threshold uj = 3 m p it gains some weight again as shown in Fig. [|b; the reader should notice, 
however, the tiny vertical scale of this Figure. 

In Fig. |5| we show the plasmon peak in the spectral function for three different values 
of rrith/T = 0.10, 0.11, 0.12. Since both the plasmon frequency and its width grow linearly 
with T, and according to Eq. ( |3. 13|) p(cu p ) ~ 1/(tilj p T), the height of the peak decreases 
with rising temperature like 1/T 2 for fixed g. This is seen in Fig. [5]. T 2 p(u p ) decreases 
accordingly with increasing coupling constant like 1/g 4 for small g. 

For small g the ratio of the damping rate to the plasma frequency is 

7(m p ,0) g 2 



LU p 647T 



(3.23) 



2 

From Eq. ( |3.22| ) we see that our resummation procedure is only applicable for A = ^ < 
~ 1. In this domain the ratio ( |3.23| ) is always smaller than about 1/10, i.e. the plasmon 
is a well-defined quasiparticle which is not overdamped. 
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B. Moving Plasmons, p / 



The evaluation of the imaginary part of the self-energy for p 7^ is a little more com- 
plicated because now two non-trivial angular integrations occur. For the g-integration we 
choose the z-axis in the direction of k — p and define t = cos(k — p, q). For the fc-integration 
we choose the z-axis along p and define u = cos(p, k). In this way the integrals in Eqs. Q3.4Q 
and ( |3~5l) can be expressed as 

The occurrence of two angles in the integrand renders the determination of the integration 
limits from the energy-conserving (^-functions more complicated. However, as shown in 
Appendix |A], similar methods can be applied as in the case p = 0. The t-integration in 
Eq. Q3.24Q can be performed by using 

dE r oik — p| , 

The condition that the zero of the energy 5-function lies inside integration domain leads to 
kinematic limits for the integrations over q, k, and u. We find (see Appendix [A]) 

9 4 T 2 1 „ /T ,\ f 1 f< dxx 2 



X — Z £ [X] 



*Mx)) f(w - e(x)) In ( f }<* 2 \ + f\ - W \ ) (3.26) 

\f(e(y 2 ))f(e( yi )+e(x)-w) 



with the same dimensionless variables ( [3.18D as before and 



z = |. (3.27) 

The limits x\, yi, and y 2 are now given by Eqs. (|A6|), (|A4|), and (|A5[), respectively, divided 
by T. In the limit p — ► (z — >• 0) Eq. ( |3.26|) correctly reduces to Eq. (|3.17|) . 
For the Landau damping contribution Im<? 2 (kSP) we obtain 

1m g 2 (u, p) = (f IT - l) J d(u, x) (1 + f(e(x))) f(w + e(x)) 

X[n {f(e(y 4 ))f(e(y 3 )-e(x)-w)) ' ^ 
12 



where 



d(u, x) = 

e(E p -u) 



1 pi \ pOO 

9{E p -uj)[6(p-uj) I du + 6(uj-p) du) dx 

1 poo ^ 1 

6{uj — E p ) ^ J du J ^ tl.r I (la I d.r 

| poo p\x+z\ 







X 



x — z e(x) 



+6(uj-E p )- de(x) d|x-z 

Z J a J\x — z\ 

if p < 3m p /2 and 



(3.29) 



d(u, x) = 
6{E p -uj) 



1 pOO 

9(E P — uj) I du dx 

uj/p Jx\ 



oo pi poo 

du I dx + du I dx 



—oj/p Jx\ 



+ 6(uj - E p ) \9(p - uj 

(pu* poo pi poo 

I du I dx + du dx 
J-l Jx\ Ju* JO 

| poo p\x+z\ 

-9{uj-E p )- / de(x) / d\ 

Z J a J\x—z\ 



X 



|x — z| e(x) 



x — z 



(3.30) 



if p > 3m p /2, respectively. In deriving these expressions from the results ( |A19[ ) and ( |A20| ) 
in Appendix [A], we used for the region uj > E p the relation 

d|k " P| = -J^P (3.31) 
du |k — p| 

to rewrite the ^-integration. Since the integrand in ( |3.28|) does not depend on u, the resulting 
integrals / d|x — z| in ( |3.29| ) and (|3.30| ) can be evaluated trivially. Again Ini(7 2 (a;,p) is non- 
zero for all uj > 0, and the full propagator has a cut along the whole real uj axis, arising 
from Landau damping as for the case p = 0. One easily checks that in the zero-momentum 
limit Eq. |T28|) reduces to (|3~T9| ). 



On the plasmon energy shell uj = , /p 2 + m 2 the plasmon damping rate can be deduced 



as 



7(VP 2 + m p> P) 



9 A T 



2567T 3 ze(z) Jo 



dx 
-L 2 



L2(0 + L2 (^)' 

/ (e-o(i-ecr 

I £(W) 2 , 



(3.32) 
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where 

This on-shell result has also recently been obtained by Jeon ]12|] by a different method. 

For the numerical evaluation of the spectral function we substitute Eqs. ( |3.26| ) and ( |3.28| ) 
into Eq. fl3.3|) and then use Eq. (|3.8|) . In Fig. || we show the spectral density at fixed 
rrith/T = 0.1 for three values of the momentum p/T. As p/T increases, the width of the 
spectral function decreases as illustrated in Fig. |7|. For large p/T the width goes to zero, 
and at the same time the height of the spectral function approaches a constant. This can 
been seen explicitly in Figs. |7] and §. This means that for large p/T the plasmon looses the 
features of a collective mode and again behaves like a non-interacting free particle! At fixed, 
non-zero momentum p the dependences on the temperature T and the coupling constant g, 
on the other hand, are qualitatively the same as at p = 0. 

From the analysis of HTL resummation in gauge theories |J it is known that there the 
summation scheme can only be guaranteed to work in the soft momentum region u,p ~ 
gT < T. A partial reason for this is that the thermal loops are momentum dependent, and 
the particular momentum dependence of the "hard" thermal loop arises as the leading term 
in a high-temperature expansion of the 1-loop self energy [|]]. Correspondingly, one might 
be inclined to distrust the above results for momenta p/T ^ 1. In <j) A theory the situation is, 
however, different in that the HTL which we resum is momentum independent to order g 2 
and g 3 ; momentum dependence enters only at order g 4 , and for weak coupling A = ^ <C 1 
as in Figs. |6], |7| it is weak on the scale m p . Therefore the essential criterium for the validity 
of the resummation procedure is the smallness of the coupling constant and not that of uj/T 
or p/T . 

This argument is also supported by a numerical investigation of the sum rule Eq. (|3.14| ). 
In Table [I] we give values for the integral up to the 3-particle threshold, 

^sum = 2/ duu p(u,p) , (3.34) 
Jo 
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for different values of m^/T (or of the coupling constant A = g 2 /2A) at p = and for 
different values of p/T at fixed m t ^/T = 0.1. From the given numbers we see that in all 
cases the sum rule is nearly exhausted by the collective plasmon mode, and that the 3-body 
decay contribution at to > E* is negligible. For large coupling constant (m t h/T > 1), where 
the resummation of hard thermal loops is no longer expected to be a good approximation, we 
begin to see appreciable violations of the sum rule. On the other hand, for m t h/T = 0.1, the 
sum rule is very well satisfied even for very large momenta, p/T ~ 100. As a consequence, 
the shape of the spectral function (in particular the interplay of the p-dependence of its 
width and height at the plasmon pole shown in Figs. |7| and is completely controlled by 
this sum rule. 



IV. CONCLUSIONS 

In this paper we studied the analytical structure of the full propagator and the proper- 
ties of the plasmon in hot 4 theory at 2-loop order, using the hard thermal loop (HTL) 
resummation scheme. 

Since in 4 theory the one-loop self-energy is completely real, damping of the plasmon 
occurs only at the 2-loop level. Its physical origin is Landau damping by scattering with 
thermal scalar field quanta. This results in a cut of the full propagator along the entire real 
frequency axis. However, for small coupling constant g the spectral function which is given 
by the discontinuity across this cut shows a strong and rather sharp plasmon peak located at 
the position expected from the 1-HTL-resummation. Its width T/u p is approximately given 
by g 2 j (327r) for plasmons at rest in the heat bath, and this width decreases even further if 
the plasmon is moving relative to the heat bath. For very fast plasmons the width goes to 
zero while the height of the plasmon peak approaches to a constant: the plasmon behaves 
again like a non-interacting free particle, i.e. it looses the features of a collective mode. 
Within the range of validity of our resummation scheme, which is defined by the smallness 

2 

of the coupling constant A = |r compared to 1, the plasmon remains a well-defined, weakly- 
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damped quasiparticle excitation for all values of the temperature. The sum rule Eq. 03.14Q 
for the spectral density is very well satisfied by the 2-loop result even for momenta much 
larger than T, and it is essentially saturated by the plasmon peak. 



ACKNOWLEDGMENTS 

We are grateful to U. Wiedemann for helpful discussions, and we would like to thank the 
referee for a careful reading of the manuscript and for pointing out an error in the original ver- 
sion of the paper. This work was supported in part by the Deutsche Forschungsgemeinschaft 
(DFG), the Bundesministerium fur Bildung und Forschung (BMBF), the National Natural 
Science Foundation of China (NSFC) and the Gesellschaft fiir Schwerionenforschung (GSI). 



APPENDIX A: KINEMATIC LIMITS FOR THE COLLISION INTEGRALS 

Here we give some technical details on how to solve for the kinematic limits of the 
integral (|3.24j) resulting from the energy conserving ^-functions. We perform the calculation 
for arbitrary external momentum p; the expressions needed in Section [III A| for plasmons at 
rest are easily obtained by setting p = 0. In this limit the integrand is independent of u, 
and the ^-integration gives a trivial factor of 2. 

We start by doing the t-integration in ( |3.24| ) using ( |3.25|) and exploiting the energy 



5-function. In order to obtain a non-zero contribution, the solution E r of the 5-function 
constraint must lie in the interval 



(|k - p| - qf + m p 2 <E r < ^Qk - p| + q) 2 + m p 2 (Al) 

defined by the limits of the t integration, — 1 < t < 1. For fixed |k — p| (or fixed k and u) 
these limits on E r generate upper and lower integration limits for q. If the latter are complex, 
the corresponding values of k and u must be excluded from the exterior integrations. 

For the 3-body decay Q3.4| ) we must solve to — E^ — E q — E r = 0, with the notation of 
Eqs. Q2.12 ), (|2.13|) . This has no solutions unless u exceeds the 3-particle threshold 
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K =m P + yP 2 + 4mg . (A2) 



With the help of the auxiliary function 



A( "' |k_Pl) = (u-E k ) 2 -ik- V ) 2 



the lower and upper integration limits for q which result from (|A1|) can be expressed as 



and 



q x = max[0, q[] with q£ = - ( y/A(u, |k - p|) - |k - p| ) ( A I ) 



1 



92 = -^A(o;,|k-p|) + |k-p|J . (A5) 

For qi, q 2 to be real, A(u>, |k — p|) must be positive. The function A is shown in Fig. [3]a. 
A(u, |k — p|) has a single zero at 

puiuj 2 — p 2 — 3 m 2 ) + ijjJiijj 2 — p 2 — 3 m 2 ) 2 — A(uj 2 — p 2 u 2 )m 2 
which for p = reduces to 



J (uj 2 — m 2 )(uj 2 — 9m 2 ) 
kl = *± — , (A7) 



a p-independent double zero at 



A;* = Juj 2 - m\ , (A* 



and a singularity in between at 

pu(uj 2 — p 2 + m 2 ) + uj^J{uj 2 — p 2 + m 2 ) 2 — 4(lu 2 — p 2 u 2 )m 2 

2(u 2 — p 2 u 



^2 - o ( . .2 „2„2^ ' 



which for p = becomes simply 



uj 2 — ml 



H=-^. (A10) 



One checks that for k > k\ the energy constraint cannot be satisfied; for k\ < k < k\ the 
function A(u, |k — p|) is negative, giving rise to complex q integration limits. Hence the 
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allowed region for k is < k < There is no constraint on u. The g-integration can be 
done analytically, giving the results presented in Section (TIT] . 

For the Landau damping term ( |3.5| ) the calculation is slightly more involved. We must 
solve uj — Ek — E q — E r = with the restriction (|Al|). For fixed |k — p| the resulting lower 
and upper limits for q can now be expressed in terms of another auxiliary function 

B(w ' |k - p|) = (^p-tk-PP (A11) 



as 



q 3 = max[0, q' 3 ] with q' 3 = - (Jb(u, |k-p|) - |k - p|) , (A12) 



and 



g 4 = - Ik - p|) + |k - p|j . (A13) 

The function B(lo, k) is illustrated in Fig. [3]b. For > E p , B(lo, |k — p|) is always positive, 
so (73, g4 are real and the ^-function can be satisfied for all k G [0, oo] and u G [—1, 1]. For 
uj < E p , B(oj, |k — p|) has a positive and real zero at 



pu(uj 2 — p 2 — 3?rip) + uj^iuj 2 — p 2 — 3m 2 ) 2 — A(lu 2 — p 2 u 2 )m 
if simultaneously u < u* with 



2 

2(ur — p z u z ) 



Aoj 2 m 2 — (uj 2 — p 2 — 3mjP 2 



w* = — * ^- P — . (A15) 

2pm p 



u* itself is real only for uj > E p where 



E P = p 2 + 4m 2 , - m p . (A16) 

For E p < uj < E p , positivity of B(uj, |k — p|) thus restricts the integration range of k to 
k > k\ if u < u* , while all k are allowed for u > u*. For uj < E p there is no upper constraint 
on u, and k is restricted to the region k > k\. 

For p = 0, there is no constraint on u, and k\ reduces to 
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\J{fnp — w 2 )(9mp — uj 2 ) 



2uj 



(A17) 



We can summarize the results so far by writing 



du / dk — > 9(E P — uj) 



9{E P -u) du 



dk 



+ 9(uj - E p 

rl 



W rOO 



1 roo ' 

/ d/c 



/■oo 
1 JO 



(A18) 



The missing lower limits of the ^-integration in two of the terms will be determined now. 

Eq. ( |A14 ) shows that for uj < p, k\ has a singularity at u — —oo/p. (At u = +uj/p, k\ 
remains finite because the numerator vanishes, too.) This implies an additional restriction 
on the lower limit of the u integration domain for the two fc^-dependent terms in Eq. ( [A 181 ) 
whenever uo < p. Which of these two terms is affected depends on the magnitude of p: For 
p < 3m p /2 we find also p < E p , so in this case the second term in Eq. (|A18|) is not affected 
by this problem because it receives contributions only from uj > E p > p. For p > 3m p /2 
we get also p > E p , and then both the first and the second term in Eq. ( |A18|) are affected. 
It turns out that the region u < —uj/p is always unphysical because either the function 
B is negative or the 5-function constraint cannot be satisfied. One thus finds, using the 
dimensionless quantities ( |3.18| ), 

p<3mp/2 



1 rl \ ;>oo 

9(E p - uj) ( 9(p - uj) I du + 9(uj-p) du) I dx 



du J dk 
9(E p - uj 



+ 9(uj - E p 

/l roo 
du J dx 



-uj/p 



du / dx + du 



-i 



-i 



dx 



(A19) 



and 



du J dk 
9(E P - uj 



p>3m p /2 



1 roo 

9(E P — uj) I du I dx 
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+ 8(u - E p ) I 8(p -u){J / du dx + J duj dx 
+ 9{uo — p)[ I du I dx + du dx 

\ J — 1 J X*_ v XI* <sO 

/l roo 
duj dx. 



(A20) 
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FIGURES 




(a) (b) 

FIG. 1. Lowest order contributions to the self energy for the resummed effective Lagrangian. 
Lines with a dot denote effective propagators with masses m t h or m p , respectively. The squared 
box denotes the additional interaction due to the mass shift. 






(a) (b) (c) 

FIG. 2. Next-to-lowest-order contributions to the self energy for the resummed effective La- 
grangian. 
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A{oi,k) 



B[(o,k) 




FIG. 3. (a) The function A(uo,k) vs. k for w > 3m p (solid line) and for uo = 3m p (dot-dashed 
line); k\ and fcg indicate the two zeroes and k\ the singularity of the solid curve, (b) The function 
B(oj, k) vs. k for to < m p (solid line), for uj = m p (dot-dashed line) and for u > m p (dashed line); 
fc| indicates the zero of the solid curve. 
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3a 1 2 3 T 



FIG. 4. The two-loop spectral function for m t h/T = 0.1. The vertical line at the peak position 
in the upper diagram indicates the spectral <5-function at one-loop order. 
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FIG. 5. The two-loop spectral function for different values of m t ^/T (or coupling constant g). 
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FIG. 6. The two-loop spectral density at fixed m^/T = 0.1 for different values ofp/T. 
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FIG. 7. (a) The width of the spectral function vs. p/T for m^/T = 0.01 (dot-dashed line), 
mth/T = 0.1 (dashed line) and rn^jT = 0.5 (solid line), (b) Close-up of (a) for small p/T. 
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FIG. 8. (a) The peak value of the spectral function vs. p/T for m t h/T = 0.01 (dot-dashed line) 
and rrip/T = 0.1 (dashed line) and m t h/T = 0.5 (solid line), (b) Close-up of (a) for small p/T. 
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TABLES 

TABLE I. The values of 7 sum for different values of m t h/T at fixed p=0 (upper half) and for 
different values oip/T at fixed mth/T = 0.1 (lower half). 



mth/T 


0.05 


0.1 


0.5 


1.0 


Turn 


0.9834 


0.9852 


0.9916 


1.0669 


p/T 


0. 


10. 


20. 


80. 


Turn 


0.9852 


0.9896 


0.9903 


0.9914 
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